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ABSTRACT 
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1.  Introduction 


Mask  matching  is  a  well  known  procedure  [l]  in  which  the  detection  of 
specific  features  in  an  image  is  carried  out  by  matching  a  set  of  templates  or 
“masks”  with  the  neighborhood  of  each  pixel  in  the  image. 

_ This  paper  describes  a  set  of  3  X  3  binary  masks  that  can  be  used  to  extract 

linear  features  from  an  image.  The  masks  are  used  to  assign  various  labels  to 
each  pixel  (each  label  corresponding  to  a  particular  mask),  and  to  associate  with 
each  label  a  set  of  confidence  measures  based  on  the  homogeneity  of  the  fore¬ 
ground  and  background  in  the  mask,  and  the  difference  between  them.  The  idea 
of  this  approach  is  to  record  a  large  set  of  useful  information  at  the  pixel  level,  in 
order  to  efficiently  make  use  of  it  at  the  later  stages  of  the  linear  feature  detec¬ 
tion  process. 

In  Section  2,  we  introduce  the  set  of  masks  and  the  confidence  measures.  In 
Section  3,  the  data  structures  and  the  algorithm  used  for  the  implementation  of 
the  masks  are  discussed. 

Section  4  contains  empirical  results  obtained  by  applying  the  masks  to 
several  image  samples.  The  principal  conclusion  drawn  from  these  results  is  that 
it  is  difficult  to  decide  which  matched  masks  are  part  of  a  linear  feature  by  sim¬ 
ple  thresholding  of  the  confidence  measures.  However,  it  is  shown  in  a  compan¬ 
ion  report  [2]  that  by  using  compatibility  relations  between  masks  at  neighboring 
pixels,  good  linear  feature  extraction  performance  can  be  obtained. 
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I5gcfigir5  gives-  a  probabilistic  analysis  of  the  frequency  of  matches  and  their 


expected  '‘robustness  ior  specific  masks  and  classes  of  masks  in  “white  noise'’ 

images.  These  results  may  help  indicate  whether  or  not  a  given  image  region 

t  9 

should  be  considered  ^interesting  ,  as  regards  frequency  of  occurrence  of 
“line-like”  masks,  for  example. 


2.  The  Masks  and  the  Confidence  Measures 


2.1.  Preliminary  Remarks 

By  applying  “mask  matching”  to  an  image,  we  mean  finding  those  NxN 
image  neighborhoods  that  match  specific  masks.  The  matching  procedure  is 
described  later  in  more  detail,  but  basically,  matching  occurs  when  there  exists  a 
threshold  dividing  the  N2  pixels  into  two  groups,  called  dark  and  bright,  such 
that  the  pattern  of  black  and  white  pixels  is  exactly  the  same  as  a  given  mask. 


Our  goal  is  to  extract  “thin”  linear  features  (1-2  pixels  wide)  from  an  image. 
For  this  task,  masks  having  N  =  1  or  N  =  2  are  obviously  too  small.  On  the 
other  hand,  higher  values  of  N  result  in  much  too  large  a  set  of  masks.  N  =  5, 
for  example,  implies  a  set  of  225  =  32 M  possible  binary  masks,  so  even  if  only 
about  20%  of  these  were  selected  as  the  subset  of  linear  feature  masks,  the 
number  of  cases  would  still  be  far  beyond  practical  levels.  Even  jV  =  4  implies 
216  =  64 K  masks,  which  is  still  very  large.  We  have  used  N  =  3  to  reduce  the 
set  of  masks  to  a  reasonable  size:  29  =  512  . 


A  representation  for  a  3  X  3  neighborhood  is  shown  in  Figure  la.  Note  that 
except  for  the  center  pixel  (referred  to  as  CP  hereafter),  the  number  assigned  to 
each  pixel  is  the  angle  it  forms  with  the  CP  (measured  counterclockwise  from  the 


positive  i -axis)  divided  by  45°. 


(*)  (b) 

Figure  1.  Representation  for  a  3x3  neighborhood,  (a)  Ordinal 
numbers  of  pixels,  (b)  The  mask  with  mask_numbtr  =  273. 

Each  of  the  possible  29  =  512  binary  masks  is  assigned  a  mask  number 
according  to  the  formula 

8 

mask_number  =  6,  2' 

i=0 

where  the  value  of  6,  for  the  ith  pixel  is  either  0  (dark)  or  1  (bright).  (See  Figure 
lb,  for  example.)  Note  that  we  have  mask_number  >  256  iff  the  mask  has  a 
bright  CP.  Moreover, 

maskjnumbtr (bright  CP)  +  maskjnumbtr {dark  CP)  =  511 

holds  for  all  pairs  of  “complementary”  masks  (i.e.,  masks  having  the  same  pat¬ 
tern  but  with  dark  and  bright  exchanged). 

A  list  of  all  possible  3X3  binary  masks  (in  increasing  order  of  their  mask 
numbers)  appears  in  Appendix  A. 


2.2.  The  Linear  Feature  Masks 


We  now  introduce  a  set  of  masks  corresponding  to  various  possible 
configurations  that  could  occur  in  a  3  X  3  neighborhood  that  intersects  a  linear 
feature.  The  masks  are  grouped  into  the  following  classes: 

Edge,  Line,  Wide-Line,  Right-Angle,  Acute-Angle,  End-Point,  Corner,  L-Line, 
T-Junction,  X-Junction  ,  and  Y-Junction. 

The  dark  CP  masks  belonging  to  each  class  are  shown  in  Figure  2.  The  remaining 
dark  CP  masks  in  the  set  are  obtained  by  rotating  each  of  the  masks  shown  in 
Figure  2  by  90°,  180°  and  270°  . 
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Figure  2.  The  dark  CP  linear  feature  masks. 
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Each  mask  defines  a  partition  of  the  3X3  neighborhood  of  the  CP  into 
dark  and  bright  pixels.  If  the  CP  is  dark,  we  will  call  the  set  of  dark  pixels  the 
foreground  and  the  set  of  bright  pixels  the  background. 

Table  1  shows  the  {/,  9  -  /}  partition  types  for  each  class  of  masks,  (where  /  and 
9  -  /  denote  the  numbers  of  foreground  and  background  pixels  respectively)  and 
gives  the  number  of  masks  of  each  partition  type. 


Table  I  -  No.  of  Masks  in  each  Class  and  Partition  Type 

Type  No. 

Class 

Partition  Types 

No.  of  Masks 

Total 

1 

Edge 

{5,4},  {6,3} 

8,  8 

D 

2 

Line 

{3,6} 

12 

D 

3 

Wide-Line 

{M} 

4 

■ 

4 

Right-Angle 

{3,6} 

8 

8 

5 

Acute- Angle 

{3,6} 

8 

8 

6 

End-Point 

{2.7} 

8 

8 

7 

Corner 

{4,5},  {7,2},  {8,1} 

4,  8,  4 

16 

8 

L-Line 

{4,5} 

8 

8 

9 

T-Junction 

{4,5},  {5,4},  {7,2} 

8,  4,  4 

16 

10 

X-Junction 

{5.4} 

2 

2 

11 

Y-Junction 

{4,5} 

8 

8 

Table  1  counts  dark  or  bright  CP  masks  only;  to  count  both,  the  numbers  should 


be  doubled.  The  total  number  of  masks  is  2  •  106  =  212. 

In  a  given  3X3  pixel  neighborhood,  let  the  nine  pixel  gray  levels  be 
0,<...<z9  .  Depending  on  how  many  of  the  z’s  are  equal,  there  can  be  up  to  8 
partitions  of  the  nine  pixels  ato  k  dark  pixels  and  9  -  k  bright  pixels,  where 
0 <  A:  <9,  obtained  by  choosing  a  threshold  between  z k  and  zk  +  x  (if  they  are  not 
equal).  Thus  any  pixel  neighborhood  can  match  up  to  8  masks  (not  all  of  which, 
however,  will  be  linear  feature  masks). 

2.3.  The  Confidence  Measures 

We  shall  use  the  following  “confidence  measures”  to  characterize  the  match 
of  a  given  pixel  neighborhood  with  a  given  mask: 

(1)  Robustness ,  defined  by  the  width  of  the  thresholding  interval,  i.e.  the 
difference  between  the  minimum  gray  level  of  a  bright  pixel  and  the  max¬ 
imum  gray  level  value  of  a  dark  pixel.  Note  that  robustness  >  0  is  a  neces¬ 
sary  precondition  for  mask  matching  to  occur. 

(2)  Contrast,  defined  as  the  difference  between  the  average  gray  levels  of  the 
dark  and  bright  pixels.  [A  more  informative  measure  can  be  obtained  by 
regarding  the  background  region  as  consisting  of  connected  components;  for 
example,  in  line  masks  there  are  two  background  components.  In  this  situa¬ 
tion  we  can  consider  separately  the  difference  between  the  average  gray  level 
of  the  foreground  pixels  and  the  pixels  belonging  to  each  background  com¬ 
ponent.] 


(3)  Fisher  distance  (see  [l],  Section  10.2.3),  defined  by  \nx  -  ji2|/(cr  x  +  , 

where  fix,  /i2  are  the  average  gray  levels  in  the  foreground  and  background, 
and  crj,  ct2  are  the  corresponding  standard  deviations.  In  Section  4  we 
explain  on  what  basis  this  measure  was  applied  to  all  the  matched  masks 
even  though  in  principle  it  should  have  been  used  only  for  single  background 
component  masks  (edge  type  masks,  for  instance). 

(4)  Homogeneity  of  the  foreground  and  each  of  the  background  connected  com¬ 
ponents.  “X-homogeneity”  (where  X  stands  for  any  connected  component 
of  pixels)  may  be  defined  as  1  -  uar(X)  /  i>armax(X)  where  var(X)  and 
t»armax(X)  denote  the  gray  level  variance  and  the  maximum  possible  gray 
level  variance  for  the  X-set  of  pixels,  respectively.  Assuming  the  gray  level 
in  the  image  takes  on  values  between  0  and  glvmAX  ,  it  can  easily  be  shown 
that 


varmaxP0 


|  1/4  9lv  max 

imi2-i 

1  iphi2 


!/4  Slv  max 


||X||  even 
otherwise 


where  ||X||  denotes  the  cardinality  of  the  X-set  of  pixels. 

(5)  Foreground  road  similarity,  defined  by  the  negative  of  the  distance  between 
the  average  foreground  gray  level  and  a  given  “expected  gray  level”  interval. 

The  measures  have  been  defined  here  for  single-band  images,  but  they  can 


easily  be  extended  to  multi-band  images. 


3.  Data  Structures  and  Algorithm 

For  several  reasons,  including  ease  of  modifiability  and  the  availability  of  a 
diverse  software  environment,  Lisp  was  the  programming  language  chosen  for  the 
implementation  of  the  proposed  system.  The  Lisp-based  data  structures  used  and 
the  algorithm  are  described  in  the  following  subsections. 

3.1.  Data  Types  and  Structures 

Using  the  extended  Franz  Lisp  environment  at  the  University  of  Maryland 
[3],  we  define  an  image  as  an  abstract  data  type  (or  Lisp  “flavor”  )  with  the  fol¬ 
lowing  “instance  variables”  or  characteristics: 

Width,  height,  number_of_bands,  band^descriptions,  datatype,  ranges,  row_scale, 
col_scale  and  data. 

Band_descnption  is  a  list  of  character  strings  giving  a  short  description  of  the 
contents  of  the  band,  one  for  each  of  the  number_of_bands  bands. 

Data_type  describes  the  type  of  data  in  the  image  arrays  (e.g.,  eight  bit  integer, 
list,  etc.).  All  bands  have  the  same  type. 

Ranges  is  a  list  of  lists  that  specify  the  range  of  data  in  each  band.  For  a  gray 
level  image,  each  element  of  the  ranges  list  is  a  list  of  two  fixnums  indicating  the 
extrema  of  the  corresponding  band  (i.e.,  the  min  and  max  gray  level  value).  For  a 
label  image,  each  element  is  a  list  of  the  possible  labels  occurring  in  that  band. 
Row_8cale  and  col_scale  are  the  numbers  of  meters  per  pixel  in  the  row  and 
column  dimensions,  respectively.  The  scale  factors  are  donums  that  are  negative 


if  unknown. 


Data  is  the  three  dimensional  array  (band  row  column)  of  data.  The  first  row  of 
data  is  the  uppermost  row  of  the  image. 

One  of  the  operations  that  can  be  performed  on  an  image  (i.e.,  one  of  the 
image  flavor  methods)  is  extracting  a  3  X  3  neighborhood  centered  around  a 
point  at  row,  col.:  (get_ 3  X  3_ neighborhood  row  col. )  The  value  returned  is  an 
ordered  list  ( center  01234567)  of  nine  elements,  where  each  element  is 
bound  to  a  sublist  consisting  of  the  gray  level  values  of  the  corresponding  num¬ 
bered  pixel  (See  Fig.  la),  in  the  different  bands  of  the  image. 

A  look-up  table  for  the  labels  assigned  to  all  possible  binary  masks  is  stored 
in  a  512-entry  mask_set  array,  where  each  element  is  bound  either  to  one  of  the 
mask  classes  or  to  nil,  depending  on  the  corresponding  entry  (or  maskjnumbtr). 

The  desired  output  of  the  mask  matching  procedure  is  an  image  (out_image 
instance)  of  the  same  size  as  the  input  image  in  which  each  non-border  pixel  is 
labeled.  For  every  labeled  pixel,  the  following  Lisp  structure  is  stored: 

((label  (mask_number  robustness  (minjthreshold)  contrast 

foregroundjhomogeneity  minjbackgroundjhomogeneity 
fis her _dis tone e  foreground_road_similarity) 

■■■) 

■■■) 

Note  that  several  labels  may  be  associated  with  each  labeled  pixel. 


3.2.  Algorithm  Outline 

To  better  understand  the  algorithm,  let  us  first  explain  the  mask  matching 
procedure  in  somewhat  more  detail.  As  indicated  earlier,  matching  occurs  when 
there  exists  a  threshold  dividing  the  pixels  of  the  3X3  neighborhood  into  two 
groups,  dark  and  bright,  such  that  the  pattern  of  dark  and  bright  is  the  same  as 
one  of  the  masks. 

To  find  such  thresholds,  we  sort  the  9  pixel  gray  levels  in  the  neighborhood. 
These  9  sorted  values  define  up  to  8  intervals  where  a  threshold  can  be  placed  to 
separate  the  neighborhood  into  two  non-empty  groups  of  pixels,  the  foreground 
and  the  background.  If  two  or  more  pixels  in  the  neighborhood  have  the  same 
gray  level,  there  will  be  fewer  than  8  intervals  in  which  to  place  the  threshold. 
(Robustness  must  be  strictly  positive.)  Each  threshold  may  give  rise  to  a  match, 
and  thus  each  pixel  may  receive  several  labels. 

To  allow  easy  modification  while  experimenting  with  various  data  samples, 
the  following  control  parameters  (in  addition  to  the  xnput^xmage  instance  and  the 
512-entry  mask_set  of  labels)  are  provided: 

•  Sub_tmage  to  examine,  namely  a  list  of  the  following  form:  (zO  y 0  ncols 
mows),  where  (zO,  yO)  denotes  the  upper  left  pixel  coordinates  of  the 
sub_\magc,  and  ncols,  mows  are  the  number  of  columns  and  rows.  Th^ 
sub_tmagc  passed  to  the  algorithm  should  not  include  border  pixels  of  the 
xnput_imagc.  Thus  if  the  xnput_xmagc  size  is  N  X  N,  the  sub_xmagc  and 
the  resulting  label  image  should  never  exceed  (,V  2)  X  {N  -  2). 


•  Minimum  acceptable  robustness  for  any  given  threshold.  The  default  value 
for  this  parameter  was  chosen  as  zero.  However,  higher  values  should  help 
eliminate  “weakly  matched”  masks. 

•  Expected  road  gray  level  interval. 

The  algorithm  can  now  be  outlined  as  follows: 

(0)  [Initialization] 

Set  xn*-(x0  +  ncols),  yn*— (yO  +  nr  ows),  x*—x0,  y*—y  0,  t*—0  and 

label*—  nil,  mask_info  *—  nil. 

(1)  [Read  3X3  Neighborhood] 

Read  the  gray  levels  of  the  3X3  neighborhood  centered  at  (  x  y  ). 

Set  n*—  (c  n 0  nl  n 2  n3  n 4  n5  n6  n 7). 

(2)  [Sort  the  Neighborhood  Gray  Levels] 

Sort  the  n  list  in  increasing  order. 

Set  sn*—(sn0  snl  sn 2  sn  3  sn4  sn5  sn6  sn  7  snS). 

(3)  [Threshold] 

If  snt  <sn(t+ 1),  partition  the  3X3  neighborhood  with  respect  to  the  <th 
threshold,  assigning  binary  coefficients  to  each  of  the  9  pixels  (for 
i=0,l,. ..,8  6,  -»— 0  if  glv{pi )< snt ,  otherwise  6,  <—  1). 

Otherwise,  goto  step  (6). 

(4)  [Evaluate  Mask  No.  and  Reference  Array] 

Evaluate  the  corresponding  maskjnumber. 

11 


Reference  the  mask_set  array;  if  label  (mask_number  )=m\,  goto  step  (6). 
Otherwise,  set  label*—  mask_set  [mask_number\. 

(5)  [Collect  Information] 

Compute  all  confidence  measures  associated  with  the  matched  masks, 
namely  robustness,  contrast,  foreground/background  homogeneities,  Fisher 
distance,  etc. 

Insert  the  above  information  into  the  mask_info  list  associated  with  the 
(x,  y)  element  of  the  out_imagc  array. 

(6)  [Increment  Threshold] 

Set  t*—t+l;  if  t  <8,  return  to  step  (3).  Otherwise,  set  £<—  0. 

(7)  [Scan  the  Image] 

Set  y*-y+V,  if  y<yn,  return  to  step  (1). 

Otherwise,  set  y«— yO,  x<— x  +  1.  If  x<xn,  return  to  step  (1). 

(8)  [Terminate] 

An  output  example,  generated  by  the  algorithm  for  a  10  X  10  sub_tmage,  is 
illustrated  in  the  following  figures.  Figures  3a  and  3b  show  the  bright  CP  masks 
in  the  corresponding  8X8  sub-image  and  (respectively)  the  robustness  and  con¬ 
trast  values  associated  with  each  mask.  (The  number  at  the  center  of  each  pixel 
denotes  its  gray  level  value.) 

In  [2],  criteria  for  the  consistency  of  pairs  of  masks  centered  at  adjacent  pix¬ 
els  are  defined.  Figure  3c  shows  the  consistent  pairs  of  masks  in  Figures  3a-b, 
linked  by  lines.  By  tracing  along  the  links,  we  see  that  the  image  contains  a 
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Figure  3a.  An  example  of  8  X  8  bright  CP  mask  output  with  associated 
robustness  values. 


^ +- * . ...w. j*  -•  .•  .  -v 


branching  linear  feature,  with  one  branch,  roughly  horizontal,  just  below  the 
center  of  the  image,  and  the  other,  roughly  diagonal,  extending  upward  and  left¬ 
ward  from  the  center. 

3.3.  Suppressing  “Weak”  Masks 

One  might  think  that  it  is  not  necessary  to  make  use  of  all  the  masks  that 
match  the  neighborhood  of  each  pixel,  but  only  of  the  “strongest”  ones,  say  those 
having  the  highest  robustness  or  contrast  values.  A  simple  experiment  shows, 
however,  that  this  is  not  true.  Figures  4-8  show,  respectively,  the  consistency 
links  obtained  when  we  keep,  at  each  pixel,  only  the  mask(s)  having 

(1)  Highest  robustness  value. 

(2)  Two  highest  robustness  values. 

(3)  Highest  contrast  value. 

(4)  Two  highest  contrast  values. 

(5)  Highest  robustness  or  highest  contrast. 

(In  cases  of  tied  values,  all  the  masks  involved  are  included.)  We  see  that  when 
this  is  done,  the  results  are  very  poor,  especially  for  methods  (1),  (3),  and  (5). 
Note  that  the  results  obtained  using  methods  (I),  (3),  and  (5)  are  very  similar, 
indicating  that  generally,  a  mask  having  the  highest  robustness  value  at  a  pixel 
also  has  the  highest  contrast  value,  and  vice  versa. 

Even  the  somewhat  better  results  obtained  using  methods  (2)  and  (4)  show 
that  elimination  of  masks  leads  to  loss  of  information.  (See  Figures  5  and  7.) 
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Figure  4.  Consistency  links  found  for  bright  CP  masks  having  the 
highest  robustness  values. 
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Thus,  it  is  not  desirable  to  suppress  any  masks  at  this  stage,  in  order  to  maintain 
all  available  information  at  the  pixel  level  for  use  in  the  linear  feature  extraction 
process. 

4.  Statistical  Studies 

It  is  of  interest  to  study  statistical  properties  of  the  confidence  measures  and 
other  parameters  associated  with  the  mask  matching  process.  The  purpose  of 
these  studies  is  to  investigate  to  what  extent  it  would  be  useful  to  make  use  of 
statistical  properties  in  the  linear  feature  detection  process.  The  studies  per¬ 
formed  treat  the  following  issues: 

•  Evaluation  of  the  frequency  of  occurrence  of  classes  of  matched  masks  in  a 
given  image. 

•  Computation  of  the  distribution  of  the  various  confidence  measures  associ¬ 
ated  with  the  matched  masks  (e.g.,  robustness,  contrast,  Fisher  distance, 
etc.),  and  of  the  corresponding  statistical  parameters  (i.e.,  the  maximum, 
mean,  median,  and  variance  values),  as  functions  of  the  class  of  matched 
masks. 

•  Comparison  between  the  statistical  results  obtained  for  the  dark  CP  masks 
and  those  for  the  bright  CP  masks,  as  a  function  of  the  class  of  matched 
masks. 

•  Comparison  between  confidence  measures  evaluated  for  matched  masks  at 
pixels  where  a  non-linear  line  detector  (described  in  Section  10.3.2  of  [l]  and 
referred  to  as  NLD  hereafter)  responds,  and  the  same  measures  evaluated  at 


pixels  where  the  above  non-linear  detector  does  not  respond. 

In  repeating  the  above  studies  for  various  image  instances,  we  want  to  find 
out  whether  or  not  there  exists  a  correlation  between  the  experimental  results 
obtained  and  the  presence  of  an  extensive  road  network  in  a  given  image.  For 
example,  it  would  be  of  interest  to  establish  that  the  higher  the  frequency  of 
occurrence  of  line-like  masks  in  a  certain  image  is,  the  more  likely  road  segments 
are  to  be  found  in  that  image. 

4.1.  Classification  of  Image  Types 

The  statistical  results  given  in  the  following  subsections  were  obtained  for 
the  sub-images  shown  in  Figures  9a-b  and  for  64  X  64  “white  noise”  images  gen¬ 
erated  using  [4]  (for  /i  =  127,  and  <7=5,  10,  15,  20,  25,  30,  35).  Figures  9a-b  show 
128  X  128  and  256  X  475  gray  level  images  respectively,  with  four  selected  win¬ 
dow  overlays  in  each  (denoted  by  wl,  w2,  ...,  w8.).  The  image  of  Figure  9a  (called 
“Greenbeltl”)  is  a  sub-image  of  the  image  shown  in  Figure  9b  (called  “Green- 
belt”),  which  was  digitized  from  an  aerial  photograph  of  the  Greenbelt,  Maryland 
area.  The  quadrants  of  Figure  9a  will  be  referred  to  below  as  “upper-left”, 
“upper-right”,  etc. 

In  order  to  relate  the  statistical  properties  to  the  presence  of  roads  in  an 
image,  we  have  classified  the  sub-images  into  the  following  three  categories: 

•  “suburban”  and  “highway”  images  -  e.g.,  windows  4,  7,  and  8. 

•  “homogeneous”  images  -  e.g.,  windows  1,  2,  3,  5,  and  6. 
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w7  w8 


Figure  9b.  “Greenbelt”  Image  (256  X  -475) 
with  four  selected  windows. 


“white  noise”  images. 


Figure  10  gives  the  gray  level  histograms  of  several  64  X  64  sub-images  of 
the  Greenbelt  area  images  and  of  two  white  noise  images.  Notice  the  general 
resemblance  of  the  histograms  of  images  of  the  same  class.  We  did  not  attempt, 
however,  to  develop  a  method  of  automatically  classifying  the  sub-images  by 
analyzing  their  histograms. 

4.2.  Statistical  Results 

4.2.1.  Frequency  of  Occurrence  of  Classes  of  Masks 

The  histograms  of  frequencies  of  occurrence  of  the  classes  of  bright  CP 
masks  (in  the  images  referred  to  in  Figure  10)  are  shown  in  Figure  11.  Again,  the 
distributions  obtained  for  images  of  the  same  type  are  characteristic.  Moreover, 
the  common  characteristics  of  this  distribution  for  images  of  different  types  (e.g., 
the  common  general  shape  and  the  maximal  frequency  obtained  for  Corner 
masks)  should  be  noticed.  However,  the  frequencies  of  occurrence  of  Edge  and 
Corner  masks  are  significantly  lower  in  the  homogeneous  and  white  noise  images, 
compared  with  the  suburban  images.  (See  Tables  3a-b  through  lOa-b  in  the  next 
subsection  for  the  exact  figures.) 

In  order  to  try  to  differentiate  between  the  various  image  types  in  terms 
more  closely  related  to  the  presence  of  roads  in  a  given  image,  we  apply  the  fol¬ 
lowing  test.  Let  the  set  of  Line,  Wide-Line,  and  L-Line  masks  be  called 
“line-like”  masks.  Given  1-2  pixel  wide  road  images,  it  is  of  interest  to  verify,  for 
example,  whether  the  frequency  of  occurrence  of  line-like  masks  increases  in 
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Figure  11.  Bright  CP  Mask  Histograms  for  64  X  64  Images. 

(a),  (h),  (c),  (d):  suburban:  upper-left,  upper-right,  lower-left,  window_7. 
(e),  (f)  homogeneous:  window^,  window_6. 

(g),  (h):  white  noise  with  a=  10,  <7=20. 


Table  2  -  Frequency  of  Occurrence  of  Line-Like  Masks  in  Various  Images 


Name 

Type 

Size 

%  of  Dark  CP 
Line-Like  Masks 

%  of  Bright  CP 
Line-Like  Masks 

1  Mask 
per  Pixel 

Total 

1  Mask 
per  Pixel 

Total 

window_4 

suburban 

20X20 

— 

24.0 

18.8 

26.2 

window_8 

suburban 

48X64 

■Hi 

15.1 

13.9 

18.5 

window_7 

suburban 

64X64 

12.5 

16.3 

14.4 

19.7 

upper-left 

suburban 

64X64 

13.7 

17.3 

13.0 

16.7 

upper-right 

suburban 

64X64 

13.8 

17.2 

13.9 

17.5 

lower-left 

suburban 

64X64 

12.5 

15.2 

14.1 

18.1 

lower- right 

suburban 

64X64 

12.1 

15.2 

12.5 

16.1 

window_l 

homogeneous 

16X16 

■■ 

SB 

10.7 

window_2 

homogeneous 

16X24 

9H 

11.7 

mmm 

window_3 

homogeneous 

16X24 

10.4 

mm 

mssm 

12.3 

window_5 

homogeneous 

64X64 

11.6 

13.7 

1  1 

11.8 

window_6 

homogeneous 

64X64 

12.1 

14.2 

mfim 

12  3 

H—127,  a= 5 

white  noise 

64X64 

■K 

17.2 

mm 

17.4 

/z=127,  o=10 

white  noise 

64X64 

Ba  1 

20.7 

WBM 

20.8 

/i  =  127,  o=15 

white  noise 

64X64 

1 

20.6 

17.2 

20.6 

H~127,  o=20 

white  noise 

64X64 

22.8 

18.2 

21.9 

^  =  127,  o=25 

white  noise 

64X64 

■H 

22.2 

17.8 

21.6 

/i  =  127,  o=30 

white  noise 

64X64 

17.9 

21.2 

18.7 

22.7 

/z=127,  o=35 

white  noise 

64X64 

19.1 

23.1 

19.7 

23.5 

suburban  images  relative  to  images  of  other  types.  The  results,  summarized  in 


Table  2,  reveal  that  this  frequency  in  fact  decreases  in  the  suburban  images 


examined,  compared  with  most  of  the  white  noise  images.  (A  theoretical  deriva¬ 


tion  of  the  frequency  of  occurrence  of  masks  in  white  noise  images  is  given  in  Sec¬ 


tion  5.)  Hence  it  is  impossible  to  conclude  that  a  given  image  region  is  “subur¬ 


ban”  based  on  its  frequency  of  occurrence  of  line-like  masks. 

> 

j 

i 


4.2.2.  Distribution  of  the  Confidence  Measures 

The  study  performed  next  computes  the  distributions  of  the  mask  confidence 
measures  (robustness,  contrast,  and  Fisher  distance)  and  their  statistical  proper¬ 
ties  (maximum,  mean,  variance,  and  median)  for  the  different  classes  of  masks  in 
various  types  of  images.  The  numerical  results  obtained  for  the  set  of  64  X  64 
images  used  in  the  previous  subsections  are  given  in  Tables  3a-b  through  lOa-b. 
(Tables  .or  additional  images  are  appended  to  this  report.)  The  robustness,  con¬ 
trast,  and  Fisher  distance  histograms  for  the  bright  CP  Line  masks  are  shown  in 
Figures  12-14.  Based  on  these  results,  the  following  observations  can  be  made. 

(1)  The  statistical  properties  are  most  similar  for  images  of  the  same  type.  (The 
numerical  values  obtained  for  the  white  noise  images  differ  from  one  another, 
but  their  statistical  properties  (as  functions  of  the  class  of  masks)  are  quali¬ 
tatively  similar.) 

(2)  The  mean  and  median  robustness  and  contrast  values  of  the  matched  masks 
of  given  classes  are  generally  greater  in  the  suburban  images  than  in  the 
homogeneous  images.  On  the  other  hand,  the  mean  and  median  Fisher  dis¬ 
tance  values  in  the  homogeneous  images  are  generally  slightly  greater  than 
in  the  suburban  images.  Also,  the  Fisher  distance  variance  over  all  matched 
masks  of  a  given  class  is  usually  smaller  in  the  homogeneous  images,  except 
in  degenerate  cases  where  the  Fisher  distance  is  infinite. 

(3)  The  values  obtained  for  the  statistical  properties  in  the  white  noise  images 
strongly  depend  on  a.  The  values  obtained  for  cr=5,  for  example,  are  similar 


to  the  values  for  the  homogeneous  images  (see  Tables  B-7a-b,  and  Tables 
7a-b,  8a-b),  whereas  the  results  for  <7=10-20  are  much  closer  to  those 
obtained  for  the  suburban  images  (see  Table  lOa-b,  Tables  3a-b  through  6a- 
b,  and  Figure  12).  Generally,  the  frequency  of  occurrence  and  the  mean  and 
median  values  of  robustness  and  contrast  of  given  classes  of  masks  in  white 
noise  images  are  increasing  functions  of  a.  (See  the  theoretical  proof  of  this 
in  Section  5.)  The  Fisher  distance  in  white  noise  images  seems  to  vary  little, 
with  respect  to  either  a  or  the  class  of  masks. 

(4)  The  contrast  histograms  of  the  bright  CP  Line  masks  in  the  white  noise 
images  strongly  resemble  the  gray  level  distributions  in  the  images  them¬ 
selves  (see  Figures  10,  13).  They  differ  greatly  in  shape  from  the  contrast  his¬ 
tograms  obtained  for  the  suburban  images;  the  contrast  of  the  Line  masks  in 
the  suburban  images  is  more  uniformly  distributed.  The  maximal  robustness 
and  contrast  values  for  the  suburban  images  are  usually  greater  than  in  the 
white  noise  images.  (See  also  the  robustness  and  contrast  histograms  of  the 
bright  CP  Line  masks,  shown  in  Figures  12-14.)  These  observations  might 
suggest  the  selection  of  masks  based  on  robustness/contrast  criteria,  but  the 
consequences  of  doing  this  are  unsatisfactory,  as  illustrated  in  Section  3.3. 

These  results  may  be  useful  in  distinguishing  between  suburban  and  homo¬ 
geneous  images,  but  do  not  seem  to  be  useful  in  distinguishing  between  suburban 
and  white  noise  images. 
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Table  5a  -  Statistics  for  Dark  Center  Pixel  Masks  in  Lower-left  Quadrant  of  Greenbeltl 
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Table  6a  -  Statistics  for  Dark  Center  Pixel  Masks  in  window  7  of  Greenbelt 
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Table  8a  -  Statistics  for  Dark  Center  Pixel  Masks  in  window  6  of  Greenbelt 
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Table  10a  -  Statistics  for  Dark  Center  Pixel  Masks  in  a  61X61  “White  Noise”  Image  with  /i  =  127, 
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Figure  13.  Contrast  Histograms  for  Bright  CP  Line  Masks  in  64  X  64 
Images. 


(a)>  (b),  (c),  (d):  suburban:  upper-left,  upper-right,  lower-left,  window_7. 
(e),  (f):  homogeneous:  window_5,  window_6. 

(g),  (h):  white  noise  with  <7=10,  <7=20. 
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Figure  14.  Fisher-Distance  Histograms  for  Bright  CP  Line  Masks  in 
64  X  64  Images. 

(a),  (b),  (c),  (d):  suburban:  upper-left,  upper-right,  lower-left,  window_7. 

(e),  (f):  homogeneous:  window_5,  window_6. 

(g),  (h):  white  noise  with  <7=10,  o—  20. 


4.2.3.  Statistics  of  Dark  vs.  Bright  CP  Masks 


Linear  features  in  an  image  include  both  bright  features  on  a  dark  back¬ 
ground  and  dark  features  on  a  bright  background.  It  is  of  interest  to  compare  the 
statistics  obtained  for  bright  and  dark  CP  masks,  which  should  respond  to  these 
two  types  of  features. 

The  results  presented  in  Tables  2  through  lOa-b  indicate  the  following  facts: 

(1)  There  is  no  significant  difference  between  the  frequencies  of  occurrence  of 
dark  and  bright  CP  masks,  even  though  some  of  the  images  seem  to  contain 
more  bright  linear  features. 

(2)  The  values  of  most  of  the  statistical  properties  are  approximately  the  same 
for  both  dark  and  bright  CP  masks  in  the  homogeneous  and  white  noise 
images. 

(3)  In  the  suburban  images,  the  mean  and  median  robustness  and  contrast  of 
the  bright  CP  masks  are  generally  greater  than  those  of  the  dark  CP  masks. 
Since  the  roads  in  the  Greenbelt  images  are  bright,  this  would  suggest  that 
it  might  be  appropriate  to  supress  weak  masks,  but  we  have  already  seen 
that  this  yields  poor  results.  A  more  interesting  suggestion  is  that  if  the 
robustness  and  contrast  values  of,  say,  the  bright  CP  masks  are  significantly 
greater  than  those  of  the  dark  CP  masks,  we  have  evidence  that  the  image  is 
not  homogeneous  or  white  noise,  and  may  contain  bright  linear  features. 
(This  ignores  the  possibility  that  an  image  may  contain  both  dark  linear 
features  on  a  bright  background  and  bright  features  on  a  dark  background.) 


However,  even  if  this  assertion  is  correct,  it  can  only  be  thought  of  as  an  aid 
in  labeling  an  image  region  as  a  whole,  but  not  as  a  means  of  detecting  indi¬ 
vidual  linear  features. 

4.2.4.  Characterization  of  “Line— Like”  Masks  at  Pixels  where  the  NLD 
Responds 

If  we  assume  that  the  nonlinear  line  detector  (NLD)  provides  “image  truth” 
for  the  location  of  linear  features,  it  is  of  interest  to  compare  the  statistics 
obtained  for  line-like  masks  at  pixels  where  the  NLD  responds  with  those 
obtained  at  pixels  where  it  does  not  respond. 

The  results,  given  in  Tables  11-12,  and  the  histograms,  shown  in  Figures 
15a-c  through  17a-c,  indicate  the  following  facts.  As  might  be  expected,  the  fre¬ 
quency  of  occurrence  of  bright  CP  line-like  masks  at  pixels  where  the  NLD 
responds  is  generally  much  greater  than  at  pixels  where  it  does  not  respond.  Simi¬ 
larly,  for  suburban  and  homogeneous  images,  the  mean  robustness  and  contrast 
values  at  pixels  where  the  NLD  responds  are  greater  than  the  values  at  pixels 
where  it  does  not  respond.  (This  is  not  true  for  white  noise  images.)  These  results 
show,  not  surprisingly,  that  the  line-like  masks  tend  to  be  matched  at  the  same 
pixels  at  which  the  NLD  responds. 


Table  11  -  Statistical  Results  for  Bright  CP  Line-Like  Masks 
at  Pixels  where  Nonlinear  Detector  Responds 

%  of  Masks  Mean  Robustness  Mean  Contrast  Mean  Fisher-distance 
Type  Size  vs.  Magnitude  vs.  Magnitude _ vs.  Magnitude _ vs.  Magnitude 
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Table  12  -  Statistical  Results  for  Bright  CP  Line-Like  Masks 
at  Pixels  where  Nonlinear  Detector  Does  Not  Respond 

%  of  Masks  Mean  Robustness  Mean  Contrast  Mean  Fisher-distance 
Type  Size  vs.  Magnitude  vs.  Magnitude  vs.  Magnitude _ vs.  Magnitude 
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Figure  15b.  Contrast  Histograms  for  Bright  CP  Line— Like  Masks  in 
Lower-left  Quadrant. 

(a) ,  (c),  (e):  at  pixels  where  nonlinear  detector  responds. 

(b) ,  (e),  (f):  at  pixels  where  nonlinear  detector  does  not  respond. 
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Figure  10a.  Robustness  Histograms  for  Bright  CP  Line-Like  Masks  in 
window_5. 

(a) ,  (c),  (e):  at  pixels  where  nonlinear  detector  responds. 

(b) ,  (e),  (f):  at  pixels  where  nonlinear  detector  does  not  respond. 
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Figure  10b.  Contrast  Histograms  for  Bright  CP  Line-Like  Masks  in 
window_5. 

(a) ,  (c),  (e):  at  pixels  where  nonlinear  detector  responds. 

(b) ,  (e),  (f):  at  pixels  where  nonlinear  detector  does  not  respond. 
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Figure  17a.  Robustness  Histograms  for  Bright  CP  Line-Like  Masks  in 
“White  Noise”  (<r=10)  Image. 

(a) ,  (c),  (e):  at  pixels  where  nonlinear  detector  responds. 

(b) ,  (e),  (f):  at  pixels  where  nonlinear  detector  does  not  respond. 
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Figure  17b.  Contrast  Histograms  for  Bright  CP  Line-Like  Masks  in 
“White  Noise”  (<7=  10)  Image. 

(a) ,  (c),  (e):  at  pixels  where  nonlinear  detector  responds. 

(b) ,  (e),  (f):  at  pixels  where  nonlinear  detector  does  not  respond. 
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4.3.  Summary  of  Experimental  Results 


(1)  The  statistics  tend  to  be  characteristic  of  images  of  a  given  type.  However, 
the  statistical  properties  of  the  white  noise  images  (of  appropriate  variance) 
overlap  those  of  the  suburban  images. 

(2)  The  frequencies  of  occurrence  of  classes  of  matched  masks  generally  do  not 
depend  strongly  on  the  type  of  the  image.  (Exceptions  are  the  End  and 
Corner  mask  classes.) 

(3)  Images  of  suburban  type  can  in  principle  be  distinguished  from  images  of 
homogeneous  type  by  the  values  of  parameters  such  as  the  mean  and  median 
robustness  or  contrast  of  the  matched  masks. 

(4)  Images  of  suburban  type  may  sometimes  be  distinguished  from  white  noise 
images  by  comparing  the  statistics  of  bright  CP  and  dark  CP  matched 


masks. 


5.  Probabilistic  Analysis 


5.1.  Motivation  and  Problem  Definition 

The  objective  of  the  analysis  in  this  section  is  to  compute  the  probabilities 
of  certain  events  (e.g.,  frequency  of  occurrence  of  a  given  mask,  expected  robust¬ 
ness  of  a  class  of  masks,  etc.)  for  a  “white  noise”  image  having  normally  distri¬ 
buted  gray  levels.  These  probabilities  can  be  used  to  identify  a  real  image  “win¬ 
dow”  as  “interesting”  in  the  sense  of  having  a  high  frequency  of  occurrence  of 
line-like  features. 

Let  glv ,  the  image  gray  level,  be  a  random  variable  having  a  Gaussian  distri¬ 
bution  with  mean  /i  and  standard  deviation  a.  Assuming  independent  pixels,  we 
want  to  evaluate  the  following: 

•  Probability  for  a  3  X  3  neighborhood  to  satisfy  a  given  mask  or  class  of 
masks. 

•  Expected  robustness  for  a  given  mask/class  of  masks. 

For  a  continuous  random  variable  glv,  the  normal  density  function  is  given 
by 

f(gtv)  = 

av2n 

Since  glv  takes  on  integer  values  only,  it  is  more  appropriate  to  treat  it  as  a  ran¬ 
dom  variable  of  discrete  type  [5].  Assuming  a  gray  scale  of  0  to  255,  the  density 
function  of  glv  is  thus 


i 


ws 


.  *\.k  ml*  m*t  , 


}{glv)  =  — J==r  2  e-ik-nf/203  ^giv-k) 
av2n  k=0 


where  6  denotes  the  Dirac  delta  or  impulse  function. 


5.2.  Theoretical  Derivation 


In  deriving  the  probability  of  a  3  X  3  neighborhood  satisfying  a  given  mask, 
we  first  observe,  based  on  the  assumption  that  the  pixels  are  statistically  indepen¬ 
dent,  that  this  probability  is  the  same  for  any  d  dark  pixel  mask.  Thus  we  need 
only  derive  the  probability  of  satisfying  a  {d,  9 -d}  mask,  without  reference  to  its 
specific  geometric  pattern. 

Let  Fglv(t)—Pr  {glv  <t }  (i.e.,  the  probability  of  a  pixel  gray  level  being  less 
than  or  equal  to  t  )  denote  the  distribution  function  of  the  random  variable  glv. 

t 

Obviously,  Pr  {glv  <t}  =  J]  f{glv—k),  and  Pr  {glv  <255}  =  1  . 

k  =0 

Let  t,  r  denote  the  threshold  and  robustness  of  a  given  {d,  9-d]  mask 
matched  by  some  3X3  neighborhood,  and  let  <t,r,d>  denote  this  statistical 
event.  The  probability  of  this  event  is  Pr  d  {glv  <  t  }Pr9~d  {glv  >  t  +  r  },  provided 
there  exists  at  least  one  “dark”  pixel  w’tli  glv  =  t  and  at  least  one  “bright”  pixel 
with  glv  =  t  +  r .  Thus 

Pr  {<  t,r,d  >  }  = 

jpr*  {glv  <  t  }  -  Prd  {glv  <  t  }  j  jpr9"1*  {glv  >  t  +r  }  -  Pr9~d  {glv  >  t+r  }  j 
or,  in  terms  of  the  density  and  distribution  functions: 


WVPV 


Pr  {<t,r,d>}  = 


{  ^,i/(  0-  [f'glv  {t)-f{glv  =  t )  j  d  j  J  [l -Fglv{t  +r  )+f(glv  =  t+r  )|  (t  +  r )  9" *  J 

Summation  of  this  probability  over  all  possible  values  for  t  and  r  (this  is  legiti¬ 
mate,  since  it  involves  addition  of  probabilities  of  disjoint  events)  yields 

255  255-r 

Pr{<d>}  =  S  E  Pr{<t,r,d>} 

r  >0  t=  0 

where  <d>  stands  for  the  event  that  a  3  X  3  neighborhood  matches  a  given 
{d,  9-d}  mask. 

It  should  be  noted  that  since  we  require  strictly  positive  robustness,  we  have 
(^)  Fr{<d>}  <  1  for  d  =  1,  2,  ...,  8  .  In  other  words,  the  probability  of  some 
3X3  neighborhood  matching  any  d  dark  pixel  mask  is  strictly  less  than  1. 

Also,  since  matching  different  masks  simultaneously  in  a  given  3X3  neigh¬ 
borhood  still  involves  disjoint  events  (i.e.,  <dx>  p]  <d2>  =  0,  provided  that 
mask_number  {d ,,  9-dj}  ^  mask_number  {d2,  9-d2}  in  case  dx  =  d2),  the  proba¬ 
bility  of  a  3  X  3  neighborhood  satisfying  a  given  class  of  masks  is 

Pr{<Cj>)  =  £  n,(.Ci,di)Pr{<d,>) 

1=1 

where  <  C;  >  denotes  the  event  that  the  neighborhood  matches  a  mask  which  is 
a  member  of  class  C},  and  n,(C;,  d, )  is  the  number  of  {</,,  9 -d,}  mask 
configurations  in  class  Cj  (See  Table  1.)  For  example,  the  probability  of  a  3  X  3 
neighborhood  satisfying  an  edge  type  mask  is 


Pr{<Edge  >}  =  8  Pr{<5>}  +  8  Pr{<6>} 


Similarly,  the  probabilities  for  the  rest  of  the  classes  can  be  evaluated. 

In  evaluating  the  expected  robustness  for  a  given  mask,  a  proper  normaliza¬ 
tion  factor  should  be  introduced.  This  implies  expressing  the  required  expected 
value  in  the  following  manner: 


E{r\<d>) 


25S 

E  r 


r>  0 


Pr  {<r,d  > } 


Pr  {<d  >} 


2S5 -r 

where  Pr  {<r,d  >}  =  Pr  {<t,r,d  >}  denotes  the  probability  that  the  neigh- 

<=  o 

borhood  matches  a  given  d  dark  pixel  mask  having  robustness  r . 


To  conclude  our  analysis,  we  now  extend  the  above  expected  value  to  a 
given  class  of  masks.  The  result  is  evaluated  by  the  following  weighted  sum: 

8 


E(r\<Cj>) 


E  «,(<?;,  <)Pr{<d,>}  P  ( r  |  <  d,  » 

i=i 


Pr{<Cj>) 


5.3.  Numerical  Results 

The  numerical  values  shown  in  the  following  tables  were  obtained  for 
^=127.5  and  a= 5,  10,  15,  20,  25,  30,  35.  (Since  the  random  variable  is  of 
discrete  type,  higher  values  for  a  are  impractical.)  Note  that  Pr  {<CU>}  = 
Pr{<Cji>}  and  P(r  |<CJ|>)  =  P(r|<Cy2>)  for  classes  of  masks  consisting 
of  single  (identical  or  complementary)  partition  types.  Thus,  the  results  for  these 
cases  are  given  for  a  combined  class  entry  form. 


Table  13  -  Probability  of  Matching  a  Given  {d,  Q-d  }  Mask  vs.  a 


.095724 

.021880 

.008889 

.005778 

.005778 

.008889 

.021880 

.095724 


<7=10 

<7=15 

<7=20 

<7=25 

8 

II 

b 

<7=35 

.103145 

.105738 

.107058 

.107857 

.108374 

.108526 

.024633 

.025635 

.026153 

.026470 

.026683 

.026837 

.010269 

.010783 

.011052 

.011217 

.011328 

.011409 

.006756 

.007125 

.007318 

.007437 

.007518 

.007576 

.006756 

.007125 

.007318 

.007437 

.007518 

.007576 

010269 

.010783 

.011052 

.011217 

.011328 

.011408 

.024633 

.025635 

.026153 

.026470 

.026683 

.026836 

.103145 

.105738 

.107058 

.107857 

.108393 

.108775 

Table  14  -  Probability  of  Matching  a  Given  Class  of  Masks  vs.  a 


Class 


Edge 

Line 

Wide-Line 

Right-Ang. 
Acute- Ang. 

End-Point 

Corner 

L-Line 

Y-Junction 

T-Junction 

X-Junction 


.117337 

.106668 

023113 

.071112 

.175040 

.581048 

.046225 

.156858 

.011556 


<7=10 

<7=15 

§ 

II 

b 

<7=25 

<7=30 

<7=35 

.136195 

.143266 

.146962 

.149233 

.150769 

.151877 

.123222 

.129399 

.132622 

.134600 

.135938 

136903 

.027024 

.028500 

.029274 

.029750 

.030072 

.030305 

.082148 

.086266 

.088415 

.089733 

.090625 

.091269 

.197064 

.205082 

.209227 

.211757 

.213464 

.214694 

.636667 

.656534 

.666731 

.672935 

.677106 

680093 

.054047 

.057000 

.058548 

.059500 

.060144 

.060610 

.179603 

.188041 

.192435 

.195128 

.196948 

.198259 

.013512 

.014250 

.014637 

.014875 

.015036 

015152 

Table  15  -  Expected  Robustness  of  a  Given  {d,  9 -d  }  Mask  vs.  a 


i 

■ 

<7=10 

(7=15 

£7=20 

cr=25 

(7=30 

<7=35 

i 

3.20248 

5.95157 

8.71042 

11.47166 

14.23372 

16.98960 

19.68012 

2 

2.28344 

4.06157 

5.85552 

7.65334 

9.45269 

11.25332 

13.06100 

3 

1.98847 

3.44697 

4.92478 

6.40733 

7.89176 

9.37728 

10.86542 

4 

1.88220 

3.22365 

4.58604 

5.95359 

7.32318 

8.69383 

10.06560 

5 

1.88220 

3.22365 

4.58604 

5.95359 

7.32318 

8.69375 

10.06437 

6 

1.98848 

3.44697 

4.92478 

6.40733 

7.89176 

9.37697 

10.86064 

7 

2.28345 

4.06158 

5.85553 

7.65335 

9.45270 

11.25231 

13.04562 

8 

3.20247 

5.95157 

8.71042 

11.47166 

14.23372 

16.98970 

19.68163 

Table  16  - 

Expected  Robustness  of  a  Given  Class  of  Ma: 

sks  VS.  <7 

Class 

Bl 

<7=10 

<7=15 

8 

II 

b 

II 

b 

o 

CO 

II 

b 

<7=35 

Edge 

Line 

1.94661 

3.35835 

4.79001 

6.22657 

7.66507 

9.10422 

10.54287 

Right-Ang. 
Acute- Ang. 

Wide-Line 

1.98847 

3.44697 

4.92478 

6.40733 

7.89176 

9.37728 

10.86542 

L-Line 

X-Junction 

Y-Junction 

1.88220 

3.22365 

4.58604 

5.95359 

7.32318 

8.69375 

10  06437 

End-Point 

2.28344 

4.06157 

5.85552 

7.65334 

9.45269 

11  25332 

13.06100 

Corner 

2.87310 

5.25078 

7.63959 

10.03116 

12.42373 

14.81250 

17.15832 

T-Junction 

2.10608 

3.68335 

5.27831 

6.87763 

8.47868 

10.08033 

11.67890 

Theoretically,  the  probability  and  expected  robustness  results  for  a  given 
{ d ,  9 -d)  mask  should  be  identical  to  those  for  a  given  {9-d,  d }  mask.  (Practi¬ 
cally,  the  results  are  identical  within  minor  inaccuracies  for  the  higher  a's  exam¬ 
ined.) 

Also,  the  agreement  between  the  theoretical  results  and  those  obtained  in 
Section  4  for  “white  noise”  images  should  be  noted.  For  example,  for  <7=5, 
Pr{<  Corner  >}  —  0.581048,  whereas  the  relative  average  frequency  of 
occurrence  of  Corner  type  masks,  in  the  corresponding  white-noise  image,  is 
approximately  59%  .  Similarly,  for  <7=20,  E(r\<Line  >)  =  6.41,  whereas  the 
average  mean  robustness  computed  for  Line  type  masks  in  the  corresponding 
white-noise  image  is  approximately  6.35  . 

Finally,  the  a  dependence  of  Pr  {<</>}  (i.e.,  {^)  Pr  {<d>}  — ►  1  as  a 

approaches  its  maximum  possible  value)  is  qualitatively  explained  by  the  follow¬ 
ing  argument.  The  higher  a,  the  less  probable  it  is  that  two  pixels  have  the  same 
gray  level  value.  Thus,  the  frequency  of  occurrence  of  masks  having  zero  robust¬ 
ness  decreases,  and  consequently,  the  probability  of  a  3  X  3  neighborhood  match¬ 
ing  any  non-zero  robustness  mask  approaches  1. 


6.  Concluding  Remarks 


In  this  paper,  we  have  introduced  a  set  of  3  X  3  binary  masks  that  can  be 
used  to  extract  linear  features  from  an  image.  By  matching  the  masks  with  the 
image,  we  assign  labels  to  each  pixel,  each  label  corresponding  to  a  particular 
mask.  We  also  compute  confidence  measures  associated  with  each  label,  based  on 
the  homogeneity  of  the  foreground  and  background  in  each  mask,  and  the 
difference  between  them.  Several  statistical  studies,  attempting  to  distinguish 
suburban  images  from  homogeneous  and  white  noise  images,  were  carried  out. 

It  was  shown  that  suppression  of  matched  masks  by  thresholding  of  the 
confidence  measures  yields  unsatisfactory  results.  It  is  preferable  to  record  all  the 
information  obtained  from  the  masks,  in  order  to  make  use  of  it  at  the  later 
stages  of  the  linear  feature  detection  process. 
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Appendix  A  -  512  3X3  Binary  Masks 
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Appendix  B  —  Statistics  of  Additional  Images. 
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Table  B  9a  -  Statistics  for  Dark  Center  Pixel  Masks  in  a  64X64  “White  Noise”  Ir  age  with  /i= 1 27,  cr=25 

Robustness  Contrast  Fisher  distance 
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Table  B-lla  -  Statistics  for  Dark  Center  Pixel  Masks  in  a  64X61  “White  Noise”  Image  with  n=l27,  cr=35 

Robustness  Contrast  Fisher  distance 

Class  #  of  Masks  - - - - - - - - - - 

Max  Mean  Var  Med  Max  Mean  Var  Med  Max  Mean  Var  Med 

;e  647  56.0  10.1  72.5  7.0  101.7  51.1  198.9  50.2  7.5  2.8  0.6  2.7 
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